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Abstract. The study by Glauber of the time-dependent statistics of the Ising 
chain is extended to the case where each spin is influenced unequally by its nearest 
neighbours. The asymmetry of the dynamics implies the failure of the detailed 
balance condition. The functional form of the rate at which an individual spin 
changes its state is constrained by the global balance condition with respect to the 
equilibrium measure of the Ising chain. The local magnetization, the equal-time 
and two-time correlation functions and the linear response to an external magnetic 
field obey linear equations which are solved explicitly. The behaviour of these 
quantities and the relation between the correlation and response functions are 
analyzed both in the stationary state and in the zero-temperature scaling regime. 
In the stationary state, a transition between two behaviours of the correlation 
function occurs when the amplitude of the asymmetry crosses a critical value, with 
the consequence that the limit fluctuation-dissipation ratio decays continuously 
from the value 1, for the equilibrium state in the absence of asymmetry, to for 
this critical value. At zero temperature, under asymmetric dynamics, the system 
loses its critical character, yet keeping many of the characteristic features of a 
coarsening system. 
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1. Introduction 

In his classical work on the time-dependent statistics of the Ising model [1] Glauber 
addresses the question of how to define the dynamics of a one-dimensional chain of 
N spins, (T„(t) = ±1, evolving from an arbitrary initial condition, in such a way that 
equilibrium is approached at large times. The spins, in contact with a heat reservoir, 
make transitions randomly between the two values ±1. The rate, or probability per 
unit time, at which an individual spin flips is also assumed to depend on the values 
of its neighbouring spins. The choice made in [1] for the rate at which the n-th spin 
flips from the value cr„ to — (T„, the other remaining fixed, is 

w{crn) = - ^70-n(o-n-i + (1-1) 

The parameter a appearing in this expression gives a time scale, while 7 is, for the 
time being, a free parameter. 

There are three possible values taken by the rate function (1.1) associated to the 
eight configurations of the group of spins {c7„; cr„_i, c7„+i} involved in the flipping of 
the spin a-„, 

^«(l-7), \(^, ^"(1 + 7)- (1-2) 

Assuming that 7 is positive, the conflgurations {+;++} and {— ; }, where fT„ 

is parallel to its neighbours, arc longer-lived than the configurations {+; } and 

{— ;-|— I-} where it is antiparallel since these configurations correspond respectively 
to the rates |a(l — 7) and \a{l + 7). The intermediate rate \a corresponds to 
the cases where the neighbouring spins arc antiparallel, {+; H — }, {+; — h}, {— ; H — }, 
{— ; — h}. The parameter 7 thus describes the tendency of spins toward ferromagnetic 
alignment and determines the equilibrium state towards which the model defined 
by the expression (1.1) of the rate relaxes. The choice 7 < would describe the 
antiferromagnetic case. 

The correspondence of the stochastic model defined so far with the Ising chain 
with energy function (or Hamiltonian) 

E{C) = -jY,^r.'Tn+U (1-3) 

n 

J denoting the coupling constant (J > 0), and C = {a\, . . . , cxAr} a spin configuration, 
is made by requesting that the rate (1.1) satisfy the detailed balance condition at 
equilibrium [1], 

P{C)w{an) = P{Cn)w{~a„), (1.4) 

where the weight P{C) is proportional to the Boltzmann factor exp(— £'(C)/T), and 
Cn denotes the configuration obtained from C by fiipping cr„. Denoting by = J/T 
the reduced coupling constant, this condition leads to the single constraint relation on 

7, 

ia(l-f7) = e^^^a(l-7), (1.5) 

which yields [1] 

7 = tanh2X. (1.6) 

With this value of the parameter 7 the rate (1.1) now ensures that the system 
prepared in an arbitrary initial condition relaxes to the equilibrium state of the Ising 
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model (1.3). The sequel of ref. [1] is devoted to the complete determination of the 
timc-dcpcndcncc of the magnetization, of the equal-time two-spin correlation function, 
of the two-time two-spin correlation function, and finally of the quantities relevant for 
a description of the model in the presence of a uniform time-varying magnetic field. 
The solvability of the model relies on the linearity of the equations giving the temporal 
evolution of the quantities of interest, itself due to the choice of the rate (1.1). 

The kinetic Ising model introduced by Glauber, defined by dynamical rules 
obeying detailed balance and relaxing from an arbitrary initial condition towards 
equilibrium, is probably the first and one of the simplest cases of a solvable model 
of nonequilibrium statistical mechanics of a strongly interacting system. This model 
has therefore, ever since, been considered as a paradigm for the investigation of new 
aspects of nonequilibrium statistical mechanics. Recent significant advances in this 
field concern phase ordering at low-temperature [2, 3], where the one-dimensional 
kinetic Ising model with symmetric dynamics at zero temperature (or more generally 
in the low temperature scaling regime) appears as the simplest example of a coarsening 
spin system [3] . Indeed, starting from a random initial condition (e.g. if the system was 
initially at equilibrium at infinite temperature), the system is unable to relax to any 
of its two ferromagnetically ordered, symmetry-related, equilibrium states, because 
the time to reach equilibrium can be arbitrarily large. On the contrary, domains of 
positive and negative magnetization grow forever, the system stays perpetually out of 
equilibrium, and, in the scaling regime, the system becomes statistically self-similar 
with only one characteristic length scale, the mean size of domains. Furthermore 
the two-time correlation and response functions only depend on the ratio of their 
arguments. The corresponding scaling functions are now well known, both for 
correlations [3, 4, 5, 6, 7, 8, 9], and for the response [9, 10, 11, 12]. Extensions to 
multispin two-time correlation and response functions can be found in [13]. 

In the present work we address the situation of a directed Ising chain where the 
flipping spin is not equally influenced by its neighbours. As for the Glauber-Ising 
model the dynamics is non conserved. The asymmetry of the dynamics means that 
the detailed balance condition is not obeyed, and therefore there is no conventional 
equilibrium -one speaks instead of a stationary state. A notable remark is that this 
stationary state does not carry any conserved macroscopic current, in contrast with 
what is usually met for the driven systems of current interest such as the Asymmetric 
Simple Exclusion Process [14], the Katz, Lebowitz and Spohn model [15], or the Zero 
Range Process [16]|. 

The definition of the model encompasses the partially asymmetric case, the totally 
asymmetric one, and the Glauber symmetric case, through a single parameter, the 
bias. For any value of the bias the stationary state of the model is the same, namely, 
the equilibrium state of the symmetric case [18]. However the relaxation properties of 
the model toward its stationary state, as well as the nature of the fluctuations in this 
stationary state, depend on the value of the bias. 

Kinetic Ising models with asymmetric conserved dynamics, such as the KLS 
model [15], represent a long-time established field of study [19]. In constrast with 
the former, very few works have been devoted to the case of Ising models with non 
conserved asymmetric dynamics. These models have been considered in the past [20] 

{ For example, in the one-dimensional KLS model [15, 17], defined as an Ising chain with conserved 
dynamics and totally asymmetric rules, only the H — bond is updated, not the — h bond. Therefore 
if the -|- spin is considered as a particle and the — spin as a hole, matter is transported in the positive 
direction. 
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and more recently [21, 18] §. Rcf. [20] gives an example of a rate function for the two- 
dimensional kinetic Ising model with totally asymmetric dynamics (up to a seemingly 
missing factor 2) leading to a Gibbsian stationary measure. However this work does 
not provide a systematic approach for the derivation of such forms of the rate, nor 
general prescription for writing its expression. Ref. [21] is concerned with the numerical 
investigation of the possible existence of a phase transition for the kinetic Ising model 
in dimensions two to five if the usual Glauber form of the rate is modified by truncating 
the local field, keeping the influential spins only. This rate neither fulfllls detailed nor 
global balance, and therefore leads to an unknown stationary measure. These two 
references motivated the work presented in [18] . The main focus was on the question: 
what is the general prescription for flnding a rate such that the stationary state be 
Gibbsian, especially for totally asymmetric dynamics? In what dimension, if any, 
does such a rate exist? The other question posed, inspired by [21], was: does the Ising 
model augmented by Glauber dynamics with the truncated rate defined above exhibit 
a phase transition to a ferromagnetic state? 

At zero temperature an alternative description of the model is in terms of diffusing 
and annihilating domain walls submitted to a bias, as explained later. This situation 
has been explored in the past [23, 24]. A discussion of the results obtained in these 
studies will be given in the course of the present work. 

The aim of the present study is to investigate, on the particular model defined 
below, the changes induced by directedness on the dynamics of the one-dimensional 
Glauber-Ising model. The agenda for the first part of the work parallels that followed 
in [1]. Wc begin by making a choice of rates for the transitions experienced by the 
flipping spin, in order to deflne the rules of the dynamics. As in the case of symmetric 
dynamics, we want to keep a correspondence between the stochastic model and the 
Ising chain with Hamiltonian (1.3). Though one can no longer rely on the detailed 
balance condition to fix the form of the rate function, it is instead possible to rely 
on the global balance condition to this end. As shown in [18], it is possible to find 
rates such that the weight of configurations at stationarity is still given by the usual 
Boltzmann-Gibbs prescription for the Hamiltonian (1.3). The model thus defined 
relaxes at long times towards a stationary state whose measure is identical to that 
of the equilibrium state of the symmetric Ising-Glaubcr model [18]. Wc shall first 
recall the results of ref. [18] concerning the determination of the rates leading to a 
stationary Boltzmann-Gibbs measure corresponding to the Hamiltonian (1.3). We 
shall then extend the study of ref. [18] to the case of a chain submitted to an external 
magnetic field. 

Having defined the process, we shall want to solve for the time dependence of 
the magnetization, and of the correlation and response fimctions for the directed Ising 
chain. Fortunately our choice of rate (2.1) still leads to linear equations, which makes 
the model solvable. The novelty comes from the introduction of another scale, through 
the occurrence of a bias (or velocity). Moreover, the absence of detailed balance, due 
to the asymmetry of the dynamics, implies the violation of the fiuctuation-dissipation 
relation, valid at equilibrium for the time-reversible kinetic Ising model with symmetric 
dynamics [1]. We shall investigate the fate of the relation between the correlation and 
response functions in such circumstances. 

§ Let us also mention a very recent publication devoted to the study of the algebraic properties of a 
disordered totally asymmetric Glauber model [22]. 



Dynamics of the directed Ising chain 



5 



Finally, paralleling the study done for the case of symmetric dynamics [9], we 
shall emphasize the low-tcmpcrature scaling behaviour of the model. For the directed 
Ising chain, coarsening is still present but it is biased. We shall present an analysis 
of this coarsening process and of the scaling behaviour of the quantities of interest 
mentioned above. While for the symmetric time-rcvcrsiblc dynamics, the fluctuation- 
dissipation relation is violated at low temperature due to the existence of an infinite 
relaxation time Teq, here this relation is already violated in the stationary state, as 
said above. We shall investigate the fluctuation-dissipation relation (in particular the 
behaviour of the fluctuation-dissipation ratio) in the zero-temperature regime where 
Teq is very large, and thus where the system is doubly a nonequilibrium system: its 
dynamics violates detailed balance and it is 'out of stationarity'. 

2. Definition of the rates for the directed Ising chain 

We consider the situation where the flipping spin fT„ is imequally influenced by its 
left and right neighbours. We shall define the rates defining the dynamical rules such 
that they lead to a stationary state measure identical to the equilibrium Boltzmann- 
Gibbs measure for the Hamiltonian (1.3), even though the dynamics is not reversible 
and therefore the detailed balance condition is violated in this stationary state. The 
general study leading to the determination of such rates is given in [18]. For the case 
of the chain with Hamiltonian (1.3) we shall merely state and explain the results. The 
case of the chain submitted to a magnetic field, which was not considered in [18], will 
be explained in greater detail. 

2.1. Rates in the absence of a magnetic field 

In the present work we use the following simple form of the rate function for the 
dynamics of the directed Ising chain, 

w{(Jn) = [1 - 70-n(P0-n-l + (1 - P)(^n+\)] , (2.1) 

where < p < 1 is a parameter which quantifies the respective influence of the 
neighbours on the flipping spin. For instance, p = (p = 1) means that ct„ only 
feels its right (left) neighbour, and p = 1/2 corresponds to the symmetric Glauber 
form (1.1). 

There are now four possible values taken by the rate function (2.1) associated to 
the eight conflgurations of the group of spins {cr„; cr„+i} involved in the flipping 
of the spin ct„, 

ia(l-7), ^a{l-jV), ^^(1 + 7^), ia(l+7), (2.2) 

where we define a velocity associated to the bias p and denoted by V = 2p — 1. 
Compared to the set of rates for symmetric dynamics given in (1.2) the only changes 
occur for the two intermediate rates where the neighbours of the flipping spin are 
antiparallel. 

The situation at zero temperature (or 7 = 1) is simpler because the rate ^a{l — j) 

corresponding to the occurrence of elementary thermal excitations (+ + H — > -\ h) 

and ( > — I — ) vanishes. Let us fix a = 1, for simplicity. The dynamics of 

the chain can entirely be described by the motion of domain walls. This is a well- 
known fact for the usual symmetric dynamics; each domain wall performs a symmetric 
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Figure 1. Space-time representation of the zero-temperature biased coarsening 
of the directed Ising chain with periodic boundary conditions (drift velocity 
V = 0.4). 



random walk, and when two domain walls encounter they annihilate. In the present 
case of directed dynamics the domain walls perform a biased random walk, with rate p 
for a step to the right and 1 — p to the left, and they again annihilate when they come 
into contact (see Figure 1). Indeed, steps to the right correspond to flipping the spin 
(T„ in configurations {-|-; — h} and {— ; H — }, with rate \{l-\-V) = p. Steps to the left 
correspond to flipping the spin C7„ in configurations {+; H — } and {— ; — h}, with rate 
^{l — V) = l—p. Ify>0 the domain walls drift away to the right, because the latter 
configurations are longer lived than the former ones. Hence, at zero temperature, 
the drift velocity of the domain walls is precisely the velocity V defined above. The 
annihilation of two domain walls corresponds to flipping the spin (t„ in conflgurations 

{+; } and {— ; ++}, with rate |a(l + 7) = 1. 

As shown below, the form of the rate (2.1) fulfills our requirement of leading to 
a stationary state measure identical to the equilibrium Boltzmann-Gibbs measure for 
the Hamiltonian (1.3). 

2.2. Constraint equations and derivation of (2.1) 

In order to understand the origin of the expression (2.1), we briefly come back on the 
main steps of the determination of the rates with the methods of [18]. The aim is 
to derive a set of constraints between the transition rates that need to be satisfied in 
order for the stationary state measure to be Gibbsian. 

The dynamics consist in flipping a spin, chosen at random, say spin n, with a rate 
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w{Cn\C), which is a more formal notation for w(fT„), corresponding to the transition 
between configurations C = {(Ti, . . . , (t„, . . . , a^} and C„ = {(Ti, . . . , — (t„, . . . , (Tat}. We 
choose periodic boundary conditions. At stationarity, the master equation expresses 
that losses are equal to gains, and reads 

P{C)Y,w{Cn\C) = Y,w{C\Cn)P{Cu), (2.3) 

n n 

where, by hypothesis, 

P(C) oc e-^(C)/T_ (2.4) 
After division of both sides by the weight P(C), oq. (2.3) can be rewritten as 

^iz;(C„|C)-«;(C|C„)e-^^/^ = 0, (2.5) 

n 

where the change in energy due to the flip reads 

AE = E{Cn) - E{C) = 2<T„ J(a„_i + <t„+i). (2.6) 

Detailed balance consists in equating each individual term appearing on both sides of 
the stationary master equation (2.5), while global balance means that this equation is 
satisfied as a whole. In other words (2.5) expresses the global balance condition on the 
rates, which are the unknown quantities of the latter. It leads to a set of constraints 
on the rates which should be satisfied in order for the stationary state to be Gibbsian 
with Hamiltonian (1.3). Before imposing any constraints the most general form of the 
rate function w((t,i) takes eight values associated to the eight possible configurations 
of the group of spins {cr„; an-iCn+i} involved in the flipping of the spin cr„, and thus 
depends a priori on eight parameters. The constraint equations as well as additional 
symmetry conditions lower this number, as ilhistratcd by what follows. 
As shown in [18], cq. (2.5) imposes the set of constraints 

W++ = e W++ 

—4,K 

w = e w 

iB+- + iu I- = w+- + w I-. (2.7) 

The compact notations Wacr' stands for ■w{an = +1;(t„_i = cr, cr„+i = a'), and Wrja' 
for w(a„ = — l;cr„_i = cr, (t„+i = a'). Note that the first two equations arc detailed 
balance conditions. These constraints lower the number of independent parameters 
to five. Requiring that the + and — spins have symmetric roles, that is imposing the 
spin symmetry condition 

Waa' = W-a-a', (2.8) 

introduces four additional constraints but reduces the set of previous constraints (2.7) 
to a single one, 

w—=e^^w++, (2.9) 

which is precisely the detailed balance condition (1.5). The most general expression 
of the rate function compatible with the constraints (2.8) and (2.9), i.e., satisfying the 
global balance condition (2.5), finally reads 

w{<Tn) = 2'^[l + cr„_icr„+i + ecr„_ia-„ 

- (7(l + <5) + e)a„a„+i], (2.10) 
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depending now on the three independent parameters, a, ^, e, and where the parameter 
7 is given by (1.6) [18]. This form of the rate corresponds to the generic situation 
where the flipping spin is unequally influenced by the left and right spins. This is the 
central result of this section. 

From this general form, three particular cases emerge. 

(i) The Glauber rate (1.1) is recovered if one imposes on (2.10) that the left and 
right spins have an equal influence on the central spin. This right-left symmetry 
therefore requires that 

w+- = «;_+, (2.11) 

a constraint which itself fixes one of the remaining unknown parameters. Carried 
in (2.10) it yields 2e + 7(1 + 5) = 0, hence 

1 7 

This expression, which now depends on two parameters, is the most general form 
of the rates for symmetric dynamics. This form, originally proposed by Glauber 
(see Appendix of ref. [1]), satisfies the detailed balance conditions (2.9) and (2.11). 
The expression (1.1) is recovered by fixing 5 = [1]. 

Note that in order to derive the form (1.1) we followed a different path than that 
chosen in [1]. We first imposed the global balance condition then symmetry 
relations, while in [1] the spin and spatial symmetry conditions are already 
encoded in (1.1), then the detailed balance condition (2.9) (i.e., (1.5)) is imposed, 
leading to (1.6), as recalled in the Introduction. 

(ii) The form of the rate (2.1) used in the present work is obtained from (2.10) with 
the particular choice of parameter 5 = 0, and using the parameterization e = —7 p. 

(iii) It is interesting to note that the form obtained from (2.1) for the totally 
asymmetric dynamics p — Q (or p = 1) is, up to the choice of time scale encoded 
in a, the unique choice at our disposal for a rate leading to a Gibbsian stationary 
state, as we now recall [18]. Assume for instance that the spin (7„ only looks to 
the right (i.e., p = 0), hence that the rates only depend on the right neighbour of 
the flipping spin: 

W++ = w I-, W-i = w (2-13) 

Carrying these conditions into (2.10), we obtain 5 = e = 0, and therefore [18] 

w{(Jn) = [1 - 7 o-„fT„+i] , (2.14) 

which is uniquely defined, up to the global time scale a. Fixing for example 
this scale by the choice a = 2cosh2iV' enables to write the rate in the compact 
exponential form 

u;(a„)=e-2^'^"^"+i. (2.15) 

The totally asymmetric dynamics is much more constrained than the partially 
asymmetric one. The existence of a rate function fulfilling the requested 
constraints (2.7), (2.8), and (2.13) was therefore not guaranteed a priori [18]. 
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2.3. Rates in the presence of a magnetic field 

In the presence of a magnetic field, corresponding to adding to the Hamiltonian (1.3) 
a term 5E = —h^^an, the methods of ref. [18] lead to the following constraint 
equations, with H = h/T, 

^++=e™«;++ 

tZ;__=e-^^+2^^i;__ 

w+- + W-+ = e'^"{w+- + W-+), (2.16) 

which are simple generalizations of (2.7). The two first equations are conditions for 
detailed balance. In the presence of a magnetic field we cannot impose the condition of 
spin symmetry (2.8). The most general expression of the rate function, corresponding 
to the general case of partially asymmetric dynamics, therefore depends on five 
parameters. This expression, that we do not write here, is the generalization of (2.10) 
for rates depending on a magnetic field. As above, imposing the additional constraints 
originating from particular choices of the symmetry of the dynamics lowers the number 
of independent parameters. We consider successively the cases of symmetric, totally 
asymmetric and finally partially asymmetric dynamics. 

(i) For symmetric dynamics, the additional constraints read 

w+_ = w I-, = w |_. (2-17) 

They lower the number of independent parameters to three and reduce the third 
equation of (2.16) to the third detailed balance condition. Instead of writing 
the associated general form for this case we give two simpler forms, both fully 
compatible with the constraint equations (2.16) and (2.17), i.e., with detailed 
balance, as follows. The first one, 

1 7 

w{cTn) = -q [1 - 2^n(i7n-l + (7„+i)] (1 - «;C7„), (2.18) 

where a fixes the time scale, and k = tanhif, was that proposed in [1]. A second 
simple expression of the rate, used in [9] , is 

w{(Tn) = [1 - C7„ tanh(ii'(£7„_i + o-„+i) + H)] . (2.19) 

(ii) For the totally asymmetric case, where for instance the flipping spin is only 
influenced by its right neigbour, the additional constraints are 

= w |_, = w , = Ou W-f = w (2.20) 

The number of independent parameters is lowered to one, the time scale. Up 
to this scale, the expression of the flipping rate is therefore once again uniquely 
determined. It reads 

W{(rn) = (1 - 70-nCT„+i) (1 - K£7„). (2.21) 

(iii) Coming back to the partially asymmetric case, we shall hereafter make the 
following choice of rate, inspired by the expressions (2.18) and (2.21), and 
interpolating between them, 

w{an) = [1 - 7(T„(pCT„_i + (1 -p)cr„+i)] (1 - Kan). (2.22) 

This form of the rate is compatible with the three constraints (2.16). In particular 
it leads to (2.21) for p = 0, as it should. 
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Let us mention that the natural generalization of (2.19), where the local field 

K{<yn-i + ""n+i) would be replaced by 2K{pan-i + (1 ^ docs not lead 

to a rate function compatible with the constraints (2.16), as long as if 7^ 0. 

Consider finally the situation where the system is subjected to a space-dependent 
arbitrary magnetic field /i„. This corresponds to adding a perturbation of the form 
SE = — hnfJn to the ferromagnetic Hamiltonian (1.3). For symmetric dynamics it 
is easily seen that extending the rate functions (2.18) and (2.19) by the replacement 
of the uniform field h by the space-dependent field /i„ is compatible with detailed 
balance. In contrast, in the presence of a space-dependent field, there is no simple 
form of the rate for asymmetric dynamics satisfying the global balance condition. This 
will not prevent the possibility of defining a response function, which only requires to 
impose an infinitesimal field to the system, as will be seen later in this work. 



3. Biased magnetization 

The first quantity we shall want to investigate is the magnetization of the chain 
Mn{t) = {an{t)). The Green function of the equations obeyed by the magnetization 
will appear as the cornerstone of much of the analysis of the next sections. 
It is an easy matter to show that M„(i) obeys the evolution equations [1] 

d(o-„) 



2{anw{an)), (3.1) 
id fixing the value of the parameter a to 1, 

= -M„(t) + 7(pM„_i(t) + (1 -p) M„+i(i)). (3.2) 



dt 

which yields, using (2.1), and fixing the value of the parameter a to 1, 
dMn{t) 



dt 

With the choice of rate (2.1), these equations are still linear, as was the case for 
symmetric dynamics [1]. They are thus easily solved by means of Laplace and Fourier 
transforms. Acting on the generic function fn{t), these transforms are defined as 
usual, as 

f^{u)= / dt/„(t)e-"*, (3.3) 
Jo 

for the temporal Laplace transform, and as 

r(9,i) = E/"(*)e"'"'' (3-4) 

n 

for the spatial Fourier transform with respect to the discrete space position n. 
Combining the two transforms we shall also make use of the double Fourier-Laplace 
transform denoted by f^^{q,u). Hereafter we shall consider an infinite system. 
Using these integral transforms, (3.2) becomes 

^^ = (^(,)-l)M-(,,i), (3.5) 

or 

uM^^{q,u) = {^{q) - l)M^^{q,u)+M^{q,t = 0), (3.6) 



yielding 



,.FT/ N M^(q,t = 0) 

^ 9. = ^ , 3.7 

u + l-(p{q) 
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where the complex function ip{q) has the alternate expressions 

(p{q) = 7 (pe~'* + (1 — p) e'^) = 7(cosg — iFsing), 

= -/\/'i--V^cos{q-qo), (tanhiq-o = I^). (3.8) 

The Green function G„(t) is the solution of (3.2) corresponding to the locally 
magnetized initial condition where the spin at the origin is pointing upwards, i.e., 
ao{t = 0) — +1, but the configuration is otherwise totaUy random. We thus have 
Gn{t = 0) = Snfi, i.e., G^{q,t = 0) = 1. Eq. (3.7) reads 

G^'^{q,u)= / (3.9) 
u+1- (p[q) 



hence 



or 



G^{q,t)=e-'^^-'^^i^\ (3.10) 



G„(t) = e-* ^ e**'(9)+™«. (3.11) 



The expression of G„(t) in Laplace space is explicit and reads 

'■^'^ dq e'"9 
27r M + 1 — 9?((j') 



G^(«) = / 
Jo 



((l + y)/(l-y))"/^ /m+1- V(u+l)2-72(l-l/2)' 



\n\ 



V(u + 1)2 -72(1-^2) \^ 7V1 - y 

(3.12) 

For symmetric Glauber dynamics with V = 0, (3.11) yields 



G'J"'{t) = e"* ^ ^ e7tcosg+m« ^ e-*7„(7i), (3.13) 

where /„ denotes the modified Bessel function of order n. A generalization of this 
result to the asymmetric case is obtained as follows. Consider the generating function 

00 

F{x,t)= J2 Gn{t)x^. (3.14) 

n=—oo 

Its explicit expression is 

F{x,t) = e-*e''*(P^+(i~P)3:), (3.15) 

as can be seen by carrying (3.14) into (3.2) and solving the resulting differential 
equation for F{x, t). Using the fact that 

^hvi^H)= J2 ^"^n(2/), (3.16) 

n=—oo 

we obtain by identification 

00 / 1 + T/\ "/2 

F{x,t)=e-' Yl ^"(1^) InhtVl^), (3.17) 

n=—oo ^ ^ 

hence the exact result 

/I I T/\ "/2 

Gn{t)=e-'(^\ /„(7iv/r^), (3.18) 
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which can easily be shown to be equivalent to (3.11), using (3.13), or alternatively 
to (3.12). The Green function is thus invariant in the simultaneous sign changes 
of n and V. Conversely, for a fixed value of V, Gn{t) is spatially asymmetric, as 
demonstrated by the expression of the time-independent ratio 

+ (3.19) 



We thus define a length scale associated to this asymmetry by 



e±VU.» = ^l±Z^ (3.20) 

where the sign ± corresponds respectively to the cases V > and y < 0. This length, 
even in V, diverges for = as ^asym ~ ±1/V and vanishes for V = ±1||. 

Coming back to the general solution of (3.2), its integral transform (3.7) yields, 
in direct space, the discrete spatial convolution 

Mnit) = M„(0) * G„(t) = J2 Mn-m{0)G^it). (3.21) 

m 

In particular, for a uniformly magnetized system, with M„(f = 0) = M, we obtain 

Mn{t) 



M 



J2Gnit)=G^{q = 0,t) = e-'/-^-, (3.22) 



at any finite temperature, where the relaxation time Teq is given by 

req=^. (3.23) 

The result obtained in (3.22) is seen to be independent of the value of the drift velocity 
V, the reason being that at zero Fourier momentum ip[q = 0) = 7 no longer depends 
on V. 

At low temperature r^q diverges exponentially fast according to Teq ~ e'^^ /2. The 
equilibrium correlation length of the Ising chain [25] , given by 

e-i/«eq =tanhX, (3.24) 

also diverges exponentially fast as ^eq ~ e^^/2, hence we have the scaling law 

Teq « 2Cq, (3.25) 

corresponding to a dynamical critical exponent z = 2, and reflecting the diffusive 
natiuc of domain growth at zero temperature, even in the presence of a non zero drift 
velocity V. 

The continuum limit estimate of Gn{t) is obtained from (3.11), by using (3.8) and 
expanding the trigonometric functions around q = 0: 

Gn{t) « e-* ^e^*(i-9'/2)+'9("-^^*), 
Jo 27r 

«^e-*(i-^)e-'^. (3.26) 
V2^ ^ ' 

II All throughout the present work the following notation is used for asymptotic equivalence; 
/(s) Ri g(x) if f{x)/g{x) 1. The sign ~ is used for asymptotic dominance by an exponential 
as in (3.31), or for indicating that two variables are of the same order of magnitude as in the example 
T ~ ^ encountered later on in the text. 
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Figure 2. Green function at zero temperature, for V = 0.4, at successive 
times t = 50,100,150,200 (from left to right). The locations of the maxima 
are at n = Vt. Continuous lines; exact result (3.18). Dashed lines: asymptotic 
estimate (3.26). 



At any positive temperature this is an exponentially damped Gaussian centered at 
n = jVt. At zero temperature this expression obeys the diffusion equation with drift 

dG„{t) _l d^Gr.it) dGr^jt) 

dt "2 9n2 ^ dn ■ ^^■^'> 
This estimate gives a faithful account of the bulk of the Green function but not of its 
tails, as depicted in Figure 2. It is valid when time t and space index n are large and 
for small velocity V, i.e., for large length scale ^asym- 

A far more accurate estimate of the Green function (3.18) is obtained by 
performing the saddle point method on the integral representation of the Bessel 
function J„(f) given in (3.13), yielding 

71^ (n2 + t2)i/4 + y;?TP j ■ (^-2^) 

This estimate is valid when time t and space n are simultaneously large. The Green 
function derived from (3.18) and (3.28) is undistinguishable from its exact value in 
Figure 2. It is maximal for n — jVt at t fixed, and, expanded around this maximum, 
yields back (3.26). From the latter estimate (3.26) one obtains the value of t where 
the Green function is maximal at n fixed and large enough: 
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if 7 = 1, and 

t « ^ ( ^ , ^ 1 , (3.30) 

in the general case. At large time the asymptotic behaviour of the Green function is 
dominated by the exponential decay 

Gn{t) ~ e-"«*, (3.31) 

with relaxation rate (or inverse characteristic time) 

ac = 1 - 7^1-^2^ (3.32) 

as can be obtained from (3.18) and (3.28), or from (3.12). We therefore confirm 
that the Gaussian approximation (3.26) is only valid for small V since it predicts 
aG = l-l(l-V^/2). 

The case of the totally asymmetric dynamics, where V = ±.1, deserves a separate 
treatment. For example, for V = 1, the Green function vanishes for n < 0, while for 
n > 0, 

This result can be obtained from the general form of the generating function (3.15), 
yielding 

F(a;,i) = e-*(^-^^). (3.34) 

It can also be obtained as the limit when ^ 1 of (3.18) using the expansion 
In{x) ~ (.T/2)"/n! for x — )• 0. Yet another means of obtaining the same result is 
by noting that, from (3.12), we have 

G^{u) = ^^V^. (3.35) 

' (u + l)"+i ^ ' 

Similarly, for V = —1, the Green function vanishes for n > 0, while for n < 0, 

^"(*) = ""(it)!'" - ^^-^^^ 
4. Equal-time correlation function 

The next quantity we shall want to study is the two-spin equal-time correlation 
function 

C{i,j,t) = {a,(t)aj{t)). (4.1) 

For a random initial condition, averaging over all the possible initial spin configurations 
and using the fact that the invariance under spatial translations along the chain is 
preserved by the dynamics, the correlation function only depends on the distance 
\j — i\ between the spins. We denote it by 

C{l,J,t)=Cn{t), (4.2) 

with n = j — i. In particular 

Coit) = 1. (4.3) 
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Again, it is easy to derive the following evolution equation [1] 

^ -2{a,a, {w{a,)+w{a,)), (4.4) 

for Cn{t), which yields a set of coupled linear differential equations of the form, 

= -2Cn{t) + l{Cn-i{t) + Cn+i{t)) (n ^ 0), (4.5) 

together with the condition (4.3), and the initial value C„(t = 0) = i5„_o- These 
equations are insensitive to the presence of the bias p as a consequence of the invariance 
of the correlation function {<7i{t)aj{t)) under the exchange of the indices i and j. 

This result is consistent with the conclusions of [23] where it was found that in 
the A + A — ^ system with driven diffusion, on which the model considered in the 
present work can be mapped at zero temperature, the local density fluctuations do 
not depend on the driving, for translationally invariant initial states. The latter are 
multispin generalizations of the equal-time correlation function studied here. 

The method given in [9] for the solution of (4.5) is now briefly recalled. We 
complete (4.5) by the corrcisponding equation for n = 0, with a time-dependent source 
v{t) in the right-hand side, to be determined in such a way that the condition (4.3) 
be fulfilled. In other words we consider the equations 

= -2Cn{t)+l{Cn-l{t)+Cn+l{t))+v{t)5n,0- (4.6) 

The solution of this equation, together with the condition (4.3), and the initial value 
Cn{i = 0) = 5nfi is given in Fourier-Laplace space by [9] 

c^\<i,u) = "^y:' = (4.7) 

^ ' u-F2-27COsg u{u + 2-2^cosq) ^ ' 



hence 



.L,„.^ _ M " + 2 - V(w + 2)2-472 



|n| 



"V^ — -) ■ (") 

We used the fact that [9] 

v^iu) = ^ V(^^ + 2)2-472- 1, (4.9) 
u 

the inverse of which is related to the correlation one site apart Ci by 

v{t) = 2{l-jCi{t)). (4.10) 

Eq. (4.8) yields limiting forms in the following two regimes of interest. 

(i) At finite temperature, for long enough times {t Teq), the correlation function 
Cn{t) converges to its equilibrium value C„,eq- We have 



CM = (4.11) 

' l-7C0sq ^ ' 

recovering thus the well-known expression Cn,eci = e-''l"l [25], where 
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(ii) At zero temperature, and more generally in the aging regime {t -C Teq), we obtain 

C„(t)«.erfc(i^), (4.13) 

where the error function erf z and the complementary error function erfc z = 
1 — erf z are defined as 



erf z = —j= \ dw e , erfc z = I du e " 

VTT Jo VT^ Jz 



(4.14) 



5. Two-time correlation function 

We now consider the two-time correlation function 

C{i,j,s,t) = {G,{s)a^{f)) , (5.1) 

where s < t is the waiting time and t = s + t the observation time. For a spatially 
homogeneous initial condition, as considered hereafter, we denote the correlation 
function by 

C{i,j,s,t) = Cr,{s,t), (5.2) 

with n = j — i. Using the same methods as above, the two-time correlation function 
is readily shown to obey the coupled linear partial differential equations 

^^^^ = -C„(s, t) + jipCn-l{s,t) + {l-p) Cn+l{s,t)), (5.3) 

for t > s, together with the initial condition 

C„(s,t = 5) =C„(s). (5.4) 

where the right-hand side is the equal-time correlation function, given in Fourier- 
Laplace space by (4.7). 

5.1. Solution of equation (5.3) 

The first argument s plays the role of a parameter in (5.3), hence these equations 
are formally identical to (3.2), with the initial condition (5.4) playing the role of the 
initial condition M„(t = 0) for the latter. The solution of (5.3), seen as an evolution 
equation in the r variable, is therefore formally identical to (3.21), and reads 

Cn{8, 8 + t) = Cn{8) * G„(r) = ^ C„_„(s)G„(r). (5.5) 

m 

In order to write down this solution more explicitly, we introduce the double 

Laplace transform of the function C„(.s, s + r), where u is conjugate to s, the transform 
being denoted by L,, and v is conjugate to r, the transform being denoted by L^: 

/•oo /"OO 

C^=L-(w,t;)= / / dsdrC„(s,s + r)e- (5.6) 
Jo Jo 

With this definition, the solution (5.5) of (5.3) reads 

v + l- ip{q} 
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i.e., using the expression (4.7) of the equal-time correlation function in Fourier-Laplace 
space, 

C^^^^^ig,u,v) = /(" + ')'-'^^' (5.8) 

' ' u{u + 2-2-/COS q){v + l- ip{q)) ^ ' 

An alternate expression is obtained as follows. Denoting by 

[u) = -J{u + 2)2 - 472 = v^' + 1, (5.9) 
u 

and noting that the Green function for symmetric dynamics, G^^™{s) is given in 
Fourier-Laplace space by 

(^sym)FL,(^^^) = ^ , (5.10) 

M + 1 — 7 COS q 

we get an expression for C„(s) in the form of a temporal convolution on s, 

C„(s) = a(s)*G^^-(2s), (5.11) 
hence from (5.5), 

Cn{s,s + T) = a{s) * ^Gn„(2s)G„(r), (5.12) 

m 

where, again, the convolution is on the variable s. 
5.2. Stationary regime 

At stationarity, letting s — oo, we find, from (5.5) or from (5.8), 

C„,stat(T) = e-"!"! * Gn{T) = ^/l^ / .(5.13) 

Jq 27r l-7cosg 

For symmetric dynamics, this expression simplifies to 

/oo 
AuG^riu), (5.14) 

which is a monotonically decreasing function. In contrast, the correlation (5.13) is not 
monotonic in r = t — s, except for n = 0, hence its derivative with respect to this time 
variable has no constant sign. 

In order to avoid this complication wc concentrate on the case n = and inspect 
the asymptotic behaviour of Go,stat(''') at large r. The two extreme cases of symmetric 
{V — 0) or totally asymmetric dynamics {V = 1) are simpler. For the former wc have 
^ofeql''') ~ as can be seen from (5.14) and (3.32). For the latter, using the 

expression (3.33) for G„(r) to compute the convolution appearing in (5.13), we obtain 
the exact expression of the correlation as 

Go,stat(r)=e-^v^. (5.15) 

The qiicstion is therefore how to interpolate between the two previous decay laws, for 
partially asymmetric dynamics (0 < y < l)Prob. The task is to compute (5.13), or 
alternatively its Laplace transform 



PrBlfe consider here the case V > Q, the extension to y < corresponding to changing the sign of V 
where appropriate. 
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The integral in this equation is, setting a = I/7 and 6 = [v + 1)/^, of the form 

dg 1 
2-K (a — cos q){h — cos q + W sin q) ' 



Z"^'' do 1 
J^K^^) = _{ „„„„ , (5-17) 



It is equal to 

6^ — a^y^ 



(5.18) 



where is a temporary compact notation for \ — V"^. The leading singularities in 
the variable b of this expression lie at bi = ^/l — V"^ and 62 = a — Vy/a^ — 1 (see the 
Appendix). They correspond to the asymptotic behaviours, at exponential order, 

Co,stat(T) ~ e-«^.=- (5.19) 

with 

ai = 1 - 7v/l - a2 = y^/l -7^- (5-20) 

The relaxation rate ai is identical to the relaxation rate ac given in (3.32). 

For \ > V > = y^l — 7^ the two singularities coexist, hence the decay is 
governed by the smaller rate, 0:2 < ai . However for V <Vc the only leading singularity 
is at 6 = 61, corresponding to the relaxation rate ai. For any positive temperature, 
there is therefore a critical velocity 14 where a transition between two behaviours 
for Constat (t) occurs. At Vc the relaxation rates are equal. These two behaviours 
interpolate between the case V = 0, where the relaxation rate is ai = 1 — 7 (unique 
singularity at 61), and the case V = 1, where the relaxation rate is a2 = ^1 — 7^ 
(unique singularity at 62), as seen previously. The same results hold if the spatial 
index n ^ 0. 

The physical interpretation of this behaviour is obtained by comparing the 
respective magnitudes of the two length scales considered earlier, namely the 
asymmetry length scale ^asym and the equilibrium correlation length ^oq- The 
condition of weak asymmetry, V < Vc, corresponds to the condition ^eq < $asym, 
the case of strong asymmetry to the condition > ^asym- 



5.3. Zero-temperature scaling regime 

General theory We can rely on the Appendix to obtain general results on the 
behaviour of the correlation function, restricting the study to Co{s,t). This quantity 
is given by (7.1). 

For symmetric dynamics, the only singularity of the integral /(a, b, 0) lies at 
6 = 61 = 1, that is 11 = 0. This corresponds to a power-law decay of the correlation at 
large time, a result which is well-known, and recovered in (5.24) below. 

For totally asymmetric dynamics, the correlation is asymptotically given by the 
product of the pure exponential e""^ by an increasing function of t, which diverges as 
e^^^, and thus sub-leading compared to the former. 

For partially asymmetric dynamics, the analysis given in the Appendix predicts 
an asymptotic exponential decay of the correlation function, at large r, with s fixed, 
with relaxation rate ai = 1 — Vl — in the whole range <V < 1. This prediction 
is confirmed by the analysis that follows. 
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Continuum limit Using the continuum limit expression of the Green function (3.26) 
we can estimate the sum appearing in (5.12) as 

G^yZ.i'^s)Gm{r) « —= / dme — e 



^/2^^^{2s + r) 

In the zero-temperature scahng regime where a^{u) ~ 2/-^, hence a{s) w 2/-^/7ri, we 
obtain, performing the convolution with a{s), 

\/2 f 1 (^-v^)" 

C„(s,s + r)w— / An , c 2(23+r-2«) . (5.21) 

^ ^ io Vw(2s + T - 2u) ^ ' 

The changes of variable, u = (2s+t)/(2 + ^;^), then z = 1/v, followed by some algebra, 
enable us to recast (5.21) into 

Q-u/2 / 



C.is, s + r).^ d. erf - j . (5.22) 

The virtue of this second form is that its derivative with respect to the time variable 
s is explicit, as will be shown below. These two equivalent expressions simplify in the 
following cases. 

(i) If T = we have 

C„(s,s)«erfc(^J^) , (5.23) 

in agreement with (4.13). 

(ii) For n = Vt, then 

2 /2s 

C„=y^(s, s + r) w — arctan \ — . (5.24) 

TT V T 

This is also the low-temperature scaling form of the two-time correlation function 
at coincidence (n = 0), in the symmetric case {V = 0) [4, 5], as can be readily 
checked in the expression (5.21), since the only dependence in n and V is through 
n — Vt. 

(iii) If r ^ s, then 

C7„(s, s + r) « erfc (^^^^^) . (5.25) 



2^ 

In particular, for n = 0, we get anomalous aging, in the sense of [26]. The two 
time variables r and s should be scaled as r ~ ^/s, which gives the bulk of the 
correlation. 

(iv) Finally if t ^ s, using the fact that the integral in (5.21) is dominated by its 
lower bound, we find 



C„(s,s + t)« J-e^ «erf(F\/i)— — . (5.26) 

V TT V\/t 

The relaxation rate of the exponential decay at large r is equal to V-^/2, which is 
the first term of the expansion of ai in V, in agreement with the prediction done 
above. 

Note that the continuum limit zero-temperature expressions (5.21) or (5.22) of 
Cn{s,t) satisfy the diffusion equation with drift (3.27). 
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In preparation of the forthcoming discussion of the fluctuation-dissipation relation, 
we need know the expressions of the derivatives of Cn{s,t) with respect to one of the 
time variables s or t, the other being kept fixed. 
Prom (5.5) we have 



dCr^is, s + r)_ dC„(s) ^ ^^(^^ ^ ^^(^^ ^ dGnir) 



ds 



ds 



ds ' 



(5.27) 



with, using (4.6), 



ds 



*Gn{T) = -2C„(s,i)+7(C„_i(s,t) + C„+i(s,i)) 
+ w(s)G„(r), 



(5.28) 



and, using (3.2), 



C„{s) * = C„(s,i) -7(pCn-i(s,t) + (1 -p)C„+i(s,i)) 

dCn{s,S + T) 



dt 



(5.29) 



Finally 



dCn{s,S + T) 

ds 



= - C„{s,t)+^{{l-p)Cn-l{s,t) +pCn+l{s,t)) 

+ t;(s)G„(T). (5.30) 

This derivative obeys the same evolution equation with respect to time t as that 
obeyed by Mn{t) and C'n{s,t). Its initial condition at t = s reads 

dCn{s,s + T) I 

dT 



\t=s 



C„(S) + j{{l-p)Cn-l{s)+pCn+l{s)) 



+ V{s)6n,0- 

We also note that 

dCn{s,S + T) dCn{s,S + T) 



(5.31) 



= v{s)Gn{T) 



ds dt 
+jV{Cn+i{s, s + t)- C„_i(s, s + r)). (5.32) 

In the zero-temperature scaling regime we have explicit expressions of these 
derivatives. Using (5.22), wc find 

dCn{s,S + T) V2 S+T 



ds 



TT {2s + T)^/ST 



e 2T- 



1 n-V(4:S + 3T) i'-vrf 
H — 7= — z^-^ — ^ e 2(2.+x) erf 



and 



/27r (2s + r)3/2 
dCr,{s,s + T) V2 Vs 



(n - Vt)^ 
Vt{2s + r) 



(5.33) 



dt 



TT (2s + T)Vr 



e 2t 



1 n + F(4s + T) (r.-v^)2 

H — = ^e 2(2s+t) erf 

V2^ (2s + r)3/2 



(n - yr)7s 
V'^(2s + r) 



(5.34) 
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0,02 




200 



Figure 3. Derivatives of Cn(s, t) with respect to s and t in the zero-temperature 
scahng regime (n = 20, s = 80, 7 = 1, V = 0.4). The derivative with respect 
to s it) is first negative (positive), then positive (negative). Continuous lines: 
asymptotic expressions (5.33) and (5.34). Dashed lines: numerical integration of 
the equations obeyed by the derivatives. 



Hence, 



dCn{s,S + T) dCn{s,S + T) 



ds 



^27r(2s + r) 



dt 



2(2a+x) erf 



72 e" 

TT 

{n - Vt)^s 
y/T{2s + T) 



(5.35) 



As depicted in Figure 3, both derivatives have no constant sign, for n fixed. The 
location where these derivatives vanish, n being fixed, can be found by the following 
analysis. Let s be large (s ^ n^r), and define the variable ^ = (n — Vt)/\/2t. Then 
from (5.33) and (5.34) we have, for C, small, 



erfC 



1 



(5.36) 



V27ryV' 

showing that the two derivatives vanish for r w n/V — l/2y^, at a fixed, large enough 
value of n. This result is similar to (3.29). 
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6. Response function and fluctuation-dissipation relation 

6.1. Overview 

Just as the two-time correlation function characterizes the spontaneous fluctuations 
of a system at thermal equilibrium, the response function characterizes how the 
system gently driven out of cqiiilibrium by a small external perturbation returns to its 
equilibrium state. In the present case of a spin system, we define the response as the 
rate of change of the limiting magnetization induced by an infinitesimal field applied 
impulsively on the system. For a time-reversible process, the fliictuation-dissipation 
theorem states that, at equilibrium, the response and correlation functions obey the 
equality 



1 dCeq(r) 

^eq(T) = -- , (6.1) 



whose physical content is that spontaneous fiuctuations and fluctuations induced by 
a small external perturbation have the same decay laws. 

Attempts to generalize this theorem to nonequilibrium situations opened several 
directions of research. 

(i) The first possible generalization consists in trying to define a temperature for a 
nonequilibrium system. A STimmary of early attempts and contributions in this 
direction is given in [27]. A natural extension of the concept of a temperature 
was then proposed in [27, 28, 29] for aging systems, such as glasses or systems 
exhibiting domain growth, as well as for systems gently driven out of equilibrium, 
through the introduction of the fluctuation-dissipation ratio X{s,t) which is the 
ratio of the response by the temporal derivative of the correlation with respect to 
the waiting time s (precise definitions are given in (6.3) and (6.31)). This quantity 
can be interpreted as giving a measure of the distance to equilibrium of the system. 
As such this definition has no further content. But for the systems of interest 
considered in [27, 28, 29] A(s, t) turns out to be a non trivial function of its two 
arguments, depending on the two times s and t through C(s, t) only. Difli'erent 
functional dependences of X{s,t) on C{s,t) correspond to difli'erent classes of out 
of equilibrium systems [27]. 

(ii) The second generalization consists in trying to obtain expressions of the response 
in terms of correlations characteristic of the non perturbed system [28, 30, 31, 32, 
33, 34, 35, 36]. When applied to non time-reversible systems, this investigation 
relies on the prior necessary question of how to define the response. 

(iii) A different viewpoint originates from the discovery that, at criticality, for systems 
undergoing critical coarsening, the limit fluctuation-dissipation ratio X^c , defined 
in (6.30), is a dimensionless amplitude ratio, and therefore a novel universal 
quantity in nonequilibrium critical dynamics [9, 37, 38]. This was demonstrated 
by the analysis of several ferromagnetic spin systems after a qiicnch from 
infinite temperature to their critical temperature, namely the one-dimensional 
Glauber-Ising chain, the two-dimensional Ising model and the spherical model 
in any dimension [9, 37, 38]. Numerous further investigations confirmed the 
universal character of the critical fiuctuation-dissipation ratio in a variety of 
systems [39, 40, 41, 42, 11]. 

In contrast, analytical and numerical studies indicate that the limit fluctuation- 
dissipation ratio vanishes throughout the low-temperature phase of coarsening 
systems [27, 43, 38]. 
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We shall be confronted with these three facets of the generalization of the 
fluctuation-dissipation theorem when examining the behaviour of the response for 
the directed Ising chain. In the present situation, the asymmetry of the dynamics 
already implies the violation of the fluctuation-dissipation theorem at stationarity. 

A preliminary question is how to define the response for such a non time-reversible 
system. As mentioned earlier, the rate function (2.22) no longer satisfies the global 
balance condition for the perturbed Hamiltonian if the applied magnetic field is space 
dependent. We shall nevertheless keep the same functional form of the rate because 
the perturbation is infinitesimal, thus the system stays infinitely close to a stationary 
state satisfying this global balance condition. 



6.2. Definition of the response function 

We define the response function in the same way as for symmetric dynamics [9]. 
We suppose that the system is subjected to a small arbitrary magnetic field /i„(i). 
We define the dynamics of the model by the rate (2.22), where k is be replaced by 
Kn{t) = tanhHn{t), with Hn{t) = hn{t)/T. The magnetization M„(t) at time t reads 

M„(t) = (a„(t)) « / duVi?„_„(«,t)/f„(u), (6.2) 

to first order in the magnetic field Hn{t), which defines the dimensionless two-time 
response function Rn-m{s,t), thus given by the functional derivative 

SM^it) 



dHm{s) 



(6.3) 

The evolution equation (3.1) yields the inhomogeneous differential equations for 

the magnetization 

= - M^{t) + j{pM^_,(t) + (l-p) M„+,{t)) 
+ K„(i)(l - -f{pC{n-l,n,t) + {l-p)Cin,n+l,t))), (6.4) 
that is, to first order in the magnetic field, 

= - M„(t) + 7(pM„_i(t) + (1 - p) M„+i(0) 

+ if„(i)(l-7Ci(i))- (6.5) 

As a consequence, the two-time response function Rn{s,t) itself obeys the linear 
differential equations 

= t) + j{pRn-i{s, t) + {l- p)Rn+i{s, t)), (6.6) 
for t > s, with the initial value 

Rnis,s) = w{s)Sn,0, (6.7) 

and with w{s) = 1 — 76*1(5) = v{s)/2 (see (4.10)), where Ci(s) is the equal-time 
correlation one site apart, given in Laplace space by (4.8). Equation (6.6), with its 
initial condition (6.7), is formally identical to (5.3) with initial condition (5.4). Hence 
its solution reads 

E„(s,s + r) = w;(s)G„(t). (6.8) 
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On view of (5.32) and (6.8) the response function can be expressed in terms of 
the correlation function C„(s, t) and its derivatives as 

1 ( dc^js.t) dCnis,t y 

^^'''^ = 2[—s or 

- ^(C„+i(s,i) - C„_i(s,i))- (6.9) 

The response is thus expressed in terms of quantities of the non perturbed system, as 
is already clear from (6.8). Generalized fluctuation-dissipation relations of this kind 
have been investigated previously [28, 32, 33, 34, 35]. The form (6.9) can be identified 
with the general expression for the response function given in [32], by noting that the 
term proportional to the drift velocity V in (6.9), which is naturally associated to the 
asymmetry due to the bias, is precisely equal to the expression named "asymmetry" 
in [32], 

Aj{s,t) = l[{c7iit)Bj{s)) - {aj{s)Bim, (6.10) 

where Bi{t) — —2aiw{(Ti) and n = j — i. 

However (6.9) does not hold in all generality. Its validity relies on a particular 
choice made for the form of the rate function. Let us take the simpler case V = 
to illustrate the point. Taking the form (2.19) for the rate, as done in [9], leads to a 
different expression for w{s), hence for As found in [9], w{s) = 1 — 7^(H-C2(s))/2. 
The corresponding response function reads 

l/aC^M) dCn{s,t) \ 7 dCi(5) 

- 2 [-^s d^J - 2 (^-^^^ 

The reason for such a discrepancy lies in the choice made in [32] for the rate in the 
presence of a fleld. It requests the latter to be of the form 

w'^ian) = M;''=°(f7„)Aw(f7„), (6.12) 

where Aw is a function of (7„, with no dependence on the neighbouring spins, which 
is not the case for the form (2.19). However, the additional term appearing in (6.11) 
vanishes identically at stationarity, and asymptotically in the zero-temperature scaling 
regime, leading in both cases to a unique definition of the response function. 



6.3. Fluctuation-dissipation relation in the presence of a uniform field 

As we noticed earlier (see (3.22)), summing upon the spatial index n i.e., taking the 
zero Fourier component of quantities, wipes out the role of the asymmetry, since 
ip{q = 0) = 7 does not depend on V. In the present case summing the two sides 
of (6.9) upon the spatial index n suppresses the contribution coming from the second 
term, proportional to V, in the right side of the equation, yielding 

pF, n .^ 1 f dC^{q = 0,s,t) dC^{q = 0,s,t) \ 

R {q = 0,s,t)=^(^ j, (6.13) 

with 

^Rn{s,t) = R^{q = 0,s,t) = ^l;(s)e-'"/^«^ (6.14) 

n 

(see (3.22)) and 

J2Cn{s,t) = C^{q = 0,s,t) = x(s)e-^/^-. (6.15) 

n 
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Figure 4. Stationary limit fluctuation-dissipation ratio versus velocity V, for 

V < Vc = — 7^ = 0.2658... (7 = tanh2). Continuous line: analytical 
prediction (6.25); dots: numerical integration of the temporal equations for Rq 
and dCo/dr. For V > Vc, Xoo.stat = 0. 



The susceptibility 

X(s)= VC„(s) = C^(<7 = 0,s), 



(6.16) 



given in Laplace space by (see (4.7)) 



X^(") = - 
u 



27 



u + 2 - 27 ^ ^ 

measures the spatial range over which ferromagnetic order has propagated at time 
s [9]. Hence C^{q — 0,s,t) and R^{q — 0,s,t) are independent of V, i.e., have the 
same forms as for symmetric dynamics. In other words, under a uniform magnetic 
field the response function is insensitive to the asymmetry of the dynamics. 



6.4- Fluctuation-dissipation relation in the stationary state 

An immediate consequence of (6.13) is that the fluctuation-dissipation theorem holds 
in the stationary state for the Fourier transforms at zero momentum, i.e., when 
applying a uniform magnetic field. 



dCLt(g = 0,T) 
dr 



(6.18) 
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Figure 5. Stationary integrated response function against stationary correlation. 
From top to bottom, V = 0, 0.2, 0.4, . . . , 1 (7 = tanh 2). 



However, in the general case, the fluctuation-dissipation theorem does not hold 
at stationarity. This can be demonstrated either using (6.9), which yields 

Rn.,tAr) = - '^^^^^f^^ - ^(a.+ l,stat(T) - C„-l.stat(T)), (6.19) 

the presence of the second term in the right side ruling out the theorem, or 
alternatively, using (6.8), yielding 

i?n,«tat (r) = G„ (r) , (6.20) 

which cannot be identified with the time derivative of C„^stat('''), given by (5.13). Note 
that, if n ^ 0, this violation is strong, because, as noticed earlier, the derivative of 
the correlation with respect to r has no constant sign, for n fixed, and therefore can 
in no way be identified to the response function, which is always positive. 

We now investigate the relative asymptotic behaviours of i?o,stat(T) and 
— dCo,stat(''')/dr. The study of the case n 7^ leads to the same conclusions. By (3.28) 
we know the asymptotic decay law of i?o,stat at large r, 

, -T(l-7\/r^) 

i?o.stat(T) « , (6.21) 

and the decay of Co,stat(''') is given by (5.19) and (5.20). 

Therefore, \iV > V^ — •\/l — 7^, the stationary limit fluctuation-dissipation ratio 
defined as 

^oo,stat = lim (6.22) 

T-!-oo -d6o,stat(T)/dT 
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vanishes, because the relaxation rate of the response ai = ac = 1 — 7\ /l — is 
always larger than the relaxation rate of the correlation a2 = V \/l — 7^. However 
for y < 14 the relaxation rates of the two quantities are both equal to a\. It turns 
out that the sub-leading corrections to the exponential decay are also the same, i.e., 
defining the amplitude Ar and Ac by 

i?o,stat (r) « AR{T)e-^(^-^^^^ (6.23) 

- ^^JWl) « Ac (r)e--(i-^^'^) , (6.24) 



dr 



we have 



Ar{t) ^. i?0,;tat(^) 



^00, Stat = lim - — 1-^ = lim — j- 



Ac'(r) "^-"it^Co^^^tatM 

1 - 7V1 - ■ 



(6.25) 



The Laplace transforms Costat(''') ^^'^ -^Ostat(^) ^'''^ respectively given by (5.16) 
and (5.18), and by (3.12), while ai is defined in (5.20). 

In stimmary, if < Vc, the stationary fluctuation-dissipation ratio has a finite 
limit at large time because the correlation and response functions have the same 
asymptotic exponential decay. For V > Vc the correlation function decays more slowly 
and Xoo,stat = 0. The unexpected result obtained here is that this effect leads to a 
transition in the behaviour of the fluctuation-dissipation ratio. As mentioned earlier, 
this transition is due to the crossing of two length scales, namely the equilibrium 
correlation length ^eq and the asymmetry length scale ^asym- The condition V < Vc 
corresponds to the condition ^eq < ^asym- The effect of the asymmetry becomes 
dominant when ^eq > Casym- 

Figure 4 gives a comparison between this prediction and the data obtained by 
a numerical integration of the temporal equations for the correlation and response 
functions. These data points were obtained by an extrapolation of a plot of the ratio 
— -Ro,stat(''')/dCo,stat(''')/dT against l/r for y < 14, or against l/\/r when approaching 
the transition. 

Another representation of the relation between correlation and response is 
obtained by using the integrated response function 

Xo(s,i) =^ duRo{u,t). (6.26) 
In the stationary state we have RQ{u,t) = i?o,stat(^ ^ "), hence 

XO,stat(T) = y du Ro,sta.t{u) = Xoo - J du Ro,sta.t{u), (6.27) 

with Xoo = -\/l — 7^7-^1 — 7^(1 — V^). In particular, at equilibrium, for V = 0, 
Xo,eq(T) = 1 - Co,eq(T). Asymptotically 

Xo,stat(T) « Xoo - M^e-«^-. (6.28) 

Hence the asymptotic ratio Xoo.stat can be extracted from the slope at the origin of a 
plot of Xo,stat(''") against Co.stat- Such a plot is depicted in figure 5, for several values 
of the velocity V and for a fixed value of the temperature. The plot was obtained by 
a numerical integration of the equations satisfied by Xn{s, t). 
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6.5. Fluctuation- dissipation ratio in the zero-temperature scaling regime 

In the presence of the two time variables s and t = s+r, we first define the asymptotic 
fluctuation-dissipation ratio [9] as 

Xas(s) = hm X„(s,s + r) = lim X„(s,t) (6.29) 

then its limit at large waiting time s, 

= lim Xas(s) = lim lim X„(s,t), (6.30) 

where the spatially dependent fluctuation-dissipation ratio reads 

General theory The study of the decay of the correlation function, done in section 5.3, 
leads to the following predictions. For symmetric dynamics the limit fluctuation- 
dissipation ratio is flnite, as is well known and as shown below. For the totally 
asymmetric dynamics we have Xas(.s) = for any waiting time s bccaiisc the decay of 
the response is faster than that of the correlation function. For partially asymmetric 
dynamics (0 < F < 1) we predict a finite asymptotic fluctuation-dissipation ratio 
Xas{s) which decays to for s — >■ oo. This last prediction is confirmed by the analysis 
of the continuum limit. 

Continuum limit In the regime of large times and V small enough such that (3.26) is 
a valid approximation of the Green function, and using the fact that w^' [u) « 1/ ^Ju 
(see (4.9)), i.e., w{s) ~ l/-\/7ri, we obtain the continuum limit of the response 



i?„(s,g + r)» \ e"'"'^'^"'', (6.32) 

7rv2sr 

which, again, is a solution of (3.27). In contrast with the two-time correlation function, 
there is no anomalous aging regime for the response because the two times s and r 
are well decoupled. 

In this regime, using (5.33) and (6.32), we find 

X„(s,s + t)-i«^^^ (6.33) 



+ - — - — —e ^(2o+^) erf 

^ (2s + r)3/2 



{n-VrWl 
^/t{2s + r) 



(6.34) 



Let us discuss some limiting cases of this expression. 

(i) For r = we obtain 

Xn{s,s)^l. (6.35) 

This result, which is known to hold in the symmetric case [9], proviso s is large, 
which is indeed the case in the continuum limit considered here, is thus extended 
to the case of a non zero asymmetry. 

(ii) More generally, if s — )• oo at fixed r, we have 
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Figure 6. Asymptotic fluctuation-dissipation ratio X^s{s) in the zero- 
temperature continuum limit against the scaling variable g = V^/s. Continuous 
line: analytical prediction (6.40), dots: numerical integration of the equations 
for correlation and response. Successively, s = 32,1/ = 0.05 = 0.282), 
s = 32,1/ = 0.1 a = 0.565), s = 16,y = 0.2 (i = 0.8), s = 32, V = 0.2 
(? = 1.13). 



hence, if n = 0, 



lim Xo(s, s + r) 



(6.37) 



where g{x) = ^/^: e^ a; erf x. 

(iii) For V and n separately equal to zero the second term in (6.33) vanishes and we 
recover the result found for symmetric dynamics [9, 10] 

^o(.,. + r)«|^. (6.38) 

The same result holds for n = Vt. For these two cases, or more generally for 
V = Q and any value of n, letting r — !■ cx) in (6.33), we obtain 

1 

2' 

i.e., = 1/2. 

(iv) For generic values of n and V ^ the expression (6.33) can be positive or negative, 
and arbitrarily large. When r — > oo, we obtain the asymptotic scaling form 

^"'^'^ " 2 + 2,g{V^) ^ 2 + erf [V ^^'^^^ 



(6.39) 
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which is a function of V^/s only, with no dependence on n. For small values of 
the scaling variable V^/s we have 

1 'W'^t 

^as(s)«-- — + •••, (6.41) 

while it converges as e~^^* to its limit value = 0. The scaling form (6.40) 
interpolates between X^o = 1/2 for V^/s = and X^^ = for V^/s oo, as 
depicted in figure 6. 

As a conclusion, the non-trivial value of the zero-temperature limit fluctuation- 
dissipation ratio Xao found for symmetric dynamics [9, 10], and characteristic of 
criticality, is lost in the presence of an asymmetry (except at the special point 
n = Vt). At zero temperature, under asymmetric dynamics, the system loses its 
critical character, yet keeping many of the characteristic features of a coarsening 
system. More precisely, the two-time quantities, namely the correlation and response 
functions, decay exponentially for asymmetric dynamics, while they have power-law 
decay for symmetric dynamics. On the other hand, the equal-time correlation function 
is insensitive to the presence of an asymmetry, and therefore all the statistics related 
to the growth of domains is unchanged. 

7. Discussion 

The model considered in the present work demonstrates how directedness, or the 
asymmetry of the dynamics, affect the propcirtics of the original Glauber- Ising model, 
which appears as a special case of the model where dynamics is symmetric. 

The two-spin equal-time correlation function turns out to be insensitive to the 
asymmetry, at any temperature. Restricting to zero temperature, where the process 
can be mapped on a one-dimensional diffusion-annihilation reaction A + A ^ with 
drift, A representing a domain wall, this result is consistent with the analysis of former 
studies [23] which predict that density fluctuations of the A-particle do not depend 
on the bias, as long as the initial state is translationally invariant. 

The effect of the asymmetry manifests itself through the emergence of a new 
scale, the asymmetry length scale ^asym- In the stationary state, two regimes for 
the behaviour of the correlation, hence of the limit fluctuation-dissipation ratio, are 
observed. These two re gimes a re separated in the velocity-temperature plane by a 
critical lino V = Vc = \/l — 7^, or equivalently ^eq = Casjm- As long as this length 
scale is greater than the equilibrium correlation length ^gq; or equivalently for values 
of the velocity less than the critical value Vc, the system departs weakly from the 
symmetric case, for which the fluctuation-dissipation theorem holds, and the response 
and correlation functions have the same asymptotic decay laws. The limit stationary 
fluctuation-dissipation ratio ^oo.stat, equal to the ratio of two amplitudes associated 
to correlation and response, is a function of temperature and velocity, taking vahies 
between 1, at V = 0, and 0, at F = Vc, and can be interpreted, in the spirit of [27], 
as giving a measure of the effective temperature of the system. Beyond the critical 
velocity Vc, i.e., when > ^asym. the effect of the asymmetry becomes dominant, and 
the correlation function has a slower decay rate than that of the response function, 
hence the limit ratio Xoo,stat vanishes. 

At zero temperature, the system is coarsening, and therefore out of stationarity. 
Yet a similar analysis predicts that, for s fixed, and for any value of V, except V = ±1, 
correlation and response have the same asymptotic decay laws at large temporal 
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separation. This can be explicitly seen in the continuum limit regime, where the 
fluctuation-dissipation ratio, in the limit r — > oo, is a function of the scaling variable 
only, which vanishes in the limit s — >■ oo. At y = ±1, corresponding to the totally 
asymmetric dynamics, the fluctuation-dissipation ratio asymptotically vanishes, when 
T — > oo, even at finite s. 

A possible extension of the present work would consist in performing the analysis 
of the behaviour of multispin two-time correlation and response functions, in order to 
investigate the role of the asymmetry, compared to what is known in the symmetric 
case [13]. At zero temperature, the two-time correlations of the energy have been 
studied in the past, in the framework of driven reaction-diffusion processes [24]. 
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Appendix 

In Laplace space the two-time correlation function Co(s, s + r) reads 

Cl^-^M = ^^±^^^,Iia,b,V) (7.1) 

where 

I{a,b,V)=r^- . (7.2) 

Jo 27r (a — cos gj(o — cos g + iK smgj 

with 

a=— — , b= . 7.3) 

27 7 

We are interested in the asymptotic behaviour of Co(s, s + t) at large r, hence in 
the singularities of the integral J(a, b, V) in the Laplace variable v, conjugate to r, or 
else in its singularities in the variable b. This integral is equal to 

1,2 2 2 

— ay 



(7.4) 



where is a temporary compact notation for 1 — V^. 
For y = it has the simpler form 

which is singular at 6 = 61 = 1. For y = 1, it has also a simpler form, reading 

J(a,6,l) = ^= = ^ ^ r^. , (7-6) 
Va^ — 1(0 — a + va — 1) 

which is singular at 6 = 62 = o — Va'^ — 1- 

In the general case, the singularities of the integral lie at 6 = 61 = \/l — for 
< y < 1, and at 6 = 62 = a - y^/a^ - 1 for K < y < 1, where K = 1/I - l/a2, 
or equivalently y = 1/a, is the value where 61 = 62. The first one corresponds to 
the vanishing of \fW^^y^, the second one to the vanishing of the quantity inside the 
brackets in the denominator of (7.4), as we now show. 

We denote this quantity by 

/(6; a, y) = ^/a^-\(b~ ay") - {a ^ b) ^b^ - y^, (7.7) 

which is a function of b, depending on the two parameters a and y. In the present 
discussion we consider a real (a > 1). Multiplying f{b;a,y) by its conjugate, 
g{b; a, y) = ^Ja^ — \ {b — ay"^) + (a — b)\/b'^ — y^, we find that the resulting expression 
vanishes for b = ±ay and for 6 = 6-1- = a ± Vy/a"^ — 1. The singularities at 6 = ±ay 
are cancelled by the numerator of I{a,b,V). As for b±, the following occurs. For a 
fixed, fib-] a,y) = for < y < 1/a, while this function is non zero for 1/a < y < 1. 
Symmetrically, g{b-;a,y) = for 1/a < y < 1, while this function is non zero for 
< y < 1/a (see Figure 7). Finally, f{b+-,a,y) is positive for < y < 1, while 
g{b^;a,y) = on the same interval. 

In summary, the only remaining singularity of I{a, b, V) coming from /(&; a, y) is 
at b-, denoted by 62 above, for 14 < y < 1. At y = 1/a, or y = 14, we have 61 = 62."'" 

+ The present discussion, done for the case V > 0, is easily extended to the case V < 0. 
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y 



Figure 7. Representation of f{b-;a,y) (in black) and g(b-.;a,y) (in red) as 
functions of y for a = 1.5. 



Coming back to the two-time correlation function Cq{s, s + r), we can conclude 
easily for the stationary case. The stationary correlation function is the limit of (7.1) 
multiplied by u, for u ^ 0, which leads to the expression (5.16) in the text, 
corresponding to a = I/7. Hence the singularity at 62 = a ~ V\/ — 1 is explicit 
and no longer depends on the Laplace variable u. 

For the case of zero temperature we proceed as follows. Asymptotically, when 
T — cx), we have, at exponential order, 

Co(s,s + T)^e^f("), (7.8) 

with g(u) = &2 — 1 = w/2 — V y^u^ /A + u. We have to invert this expression with 
respect to the Laplace variable u. Since t is large, the saddle point method can be 
applied to the inverse Laplace transform, yielding g'(u) — ~s/t — > 0. The relevant 
solution Mo = 2(l/-\/l — V"^ — 1) corresponds to the minimum of g{u). Carried in (7.8) 
this yields, at exponential order, 

Co(s,s + T)^e--(i-^^^), (7.9) 

which corresponds precisely to the singularity 61. The conclusion is that for partially 
symmetric dynamics (0 < y < 1) the asymptotic decay of the correlation function at 
large r with s fixed is governed by a single relaxation rate ai = 1 — Vl — V^. 

Finally, at finite temperature, the very same analysis leads to the same conclusion 
with ai = 1 — 7V1 — V^. 
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The case of totally asymmetric dynamics must be treated separately. We have 
from (7.1) and (7.6), with 7=1, 

cl-'^^M = r^TTi^v (^-^^^ 

U{V — u/2 + y/U''/4: + U) 

yielding 

C^'{u,T) = ^e^^^-'/?+^l (7.11) 

In this case, the minimum of g{u) = u/2 - ^u^ /4 + u is reached at infinity as can be 
seen on the expression of uq. Expanding g(u) for u — >■ oo gives g(u) « — 1 + Xju. By 
inversion of C^'{u,t') k, e~'^(^~^/")/u with respect to u we finally obtain 

Co(s, s + r) « e-^/o(2\/i7) ~ e'^+^v^. (7.12) 

Let us mention another expression of the correlation for totally asymmetric 
dynamics, 

Co(s,s + r) = C„(s)G„(t) = e-- '^^^^^^V ^^'^"^^ 

m>0 m>0 

Its Laplace transform with respect to s yields back (7.11), using (4.8). 
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